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Abstract 



We present a theoretical framework which establishes how the core radius of a star cluster varies with 
the mass of an assumed central black hole. Our result is that r c /rh oc (Mb/j/M) 3 / 4 when the system is well 
\^ | relaxed. The theory compares favourably with a number of simulations of this problem, which extend to 

black hole masses of order 10% of the cluster mass. Though strictly limited as yet to clusters with stars of 
equal mass, our conclusion strengthens the view that clusters with large core radii are the most promising 
candidates in which to find a massive black hole. 
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r-j 1. Introduction ation is mimicked by a suitably crafted form for ther- 
IL" ■ mal conductivity in a self-gravitating gas (Amaro-Seoane 
. £h ! Though the existence of intermediate-mass black holes et al., 2004); a tree code (Arabadjis, 1997); and, most 
in star clusters remains controversial, our theoretical un- recently, direct A-body simulations (Baumgardt et al., 
J-j I derstanding of the problem has advanced on two fronts. 2004a, b, 2005; Trenti et al., 2007). 
. . . i First, important conclusions have been reached from In many of these studies, emphasis is given to the de- 
purely theoretical considerations, including the density tails of the cusp and the growth of the black hole, and less 
profile of the cusp surrounding the black hole (Bahcall attention is paid to the evolution of the star cluster. In 
& Wolf, 1976; Shapiro & Lightman, 1976). This can be this letter we shall show that rather simple considerations 
understood as the response of the stellar distribution to allow us to predict also the evolution of the structure of 
the steady transport of stars into the vicinity of the black the cluster, in particular its core radius. The basic idea 
hole, where they are tidally disrupted, contributing to is a familiar one in the stellar dynamics of star clusters, 
the growth of the black hole. Two-body relaxation is where it was introduced by Hcnon (1975). This idea is 
the vital process which controls the rate at which this that the flux of energy from the centre of a star cluster 
flux can be sustained. Second, a succession of Simula- must reach an equilibrium with the flux (by relaxation) 
tions have added much detail to the general theoretical across the outer parts of the cluster, conventionally taken 
picture. These simulations have been based on a variety to be the half-mass radius. If too much energy is generated 
of techniques: Monte Carlo methods based on a Fokker- the core of the star cluster must expand to quench the gen- 
Planck treatment of relaxation with an anisotropic distri- eration of energy there, no matter what is the mechanism 
bution of velocities (Shapiro & Marchant, 1978; Duncan & of energy generation; while insufficient generation of en- 
Shapiro, 1982; Frcitag & Bcnz, 2002); finite-difference so- ergy leads to the familiar process of core collapse. For the 
lution of the Fokkcr-Planck equation for both anisotropic case of a central black hole, energy is generated as stars 
(Cohn & Kulsrud, 1978) and isotropic distribution func- fall by relaxation towards the radius around the black hole 
tions (Murphy et al., 1991); gas models, in which relax- 
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at which they are disrupted 1 . Though the link between 
the energy generated in the cusp and the expansion of 
the system has been discussed previously (see Sec. 5), we 
believe that the theoretical estimate of the core radius in 
this letter (eq.(5)) is novel. 

Here is an outline of the letter. In the following section 
we review the basic results on the flow of energy through 
the cusp around the black hole. Then we apply Henon's 
argument, which establishes the way in which the radius 
of the core varies with the mass of the black hole. Next 
we compare our prediction with existing data from simu- 
lations, and finally discuss the place of our result within 
the literature on this subject. 

2. The Energy Flux from a Central Black Hole 

We consider a cluster of mass M containing a central 
black hole of mass Mbh- The black hole is surrounded by 
a cusp, which merges at its edge into a core of nearly con- 
stant density (Fig.l). Beyond the core radius the density 
falls off again in the halo of the star cluster. 

Inside the cusp, the flux of energy at radius r is 
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where p is the stellar density, v 2 is the mean square stellar 
velocity, and t r is the relaxation time. This is given by 



t 



v 



-, where In A is the Coulomb logarithm and 



G 2 mp\n A ' 
m is the individual stellar mass. Hence 

G 2 mp 2 r 3 In A 

v 

In the part of the cluster where the potential is dominated 
by the black hole, i.e. r < ~— , where v c is the velocity 



dispersion in the core of the cluster, we have v 
and so 



GMbh 
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If the flux is independent of r (Lightman & Shapiro, 1977) 
we get p cx r~ 7 / 4 . 

At the edge of the cusp we have v 2 ~ v 2 , and so the 
radius of the cusp is 



Tcusp 



GM bh 



(1) 



At this radius we also have p~p c , where p c is the density 
in the core, and so 
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(2) 



1 Our arguments also apply, in principle, to a purely classi- 
cal, point-mass idealisation, in which stars simply accumulate 
around the black hole at small distances. 



3. The Radius of the Core 

In steady post-collapse expansion this must balance the 
energy flux at the half-mass radius, which is 



, Mv 2 h 

Oh ~ 
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(3) 
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where the subscript h denotes conditions at the half-mass 

M 

radius. Estimating ph ~ —r and equating £h to £ , we 
obtain 

M 2 G^plMl 



v h rl v\ 

If conditions are approximately isothermal between the 
core radius, r c , and the half-mass radius, we can estimate 



^4 and «g~t*~^, whence 

r h r c 
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(5) 



This scaling of core radius with black hole mass is ex- 
pected to be approximately valid provided that the result- 
ing core radius exceeds that of the cusp around the black 
hole (eq.l). Using the above estimate for v 2 , it follows 
that 



r cusp I Mbh 



1/4 



(6) 



For a sufficiently massive black hole, the density profile 
may show no sign of a core radius. This conclusion was 
already reached by Marchant & Shapiro (1979), though 
they modelled systems with a core radius independent of 
the black hole mass. A fit to numerical data (Sec. 4), how- 
ever, yields a coefficient of about 0.7 in eq.(6), and we 
find that, even up to a black hole mass of order 10% of 
the cluster mass, the radius of the cusp is less than half 
the core radius. On the other hand unless Mbh 1S large 
enough, the cusp contains too few stars to provide a sig- 
nature of the presence of the black hole (Baumgardt et 
al., 2005). 

4. Comparison with Simulations 

Baumgardt et al. (2004a) have carried out a series of 
A^-body simulations with stars of equal mass and several 
different values for the initial black hole mass. In addi- 
tion, unpublished data from some of these runs allow us 
to measure the core radius for intermediate values of the 
black hole mass, which increases during the course of each 
simulation. The measurement of core radius is not en- 
tirely straightforward, however. While the W-body code 
itself gives a current value of the "core radius" (Aarseth, 
2003), this is based on a density- weighted average of the 
distances of stars from the density centre. Because the 
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Fig. 1. Example of a fit of eq.(7) to the density profile. Data 
outside r = 2 were not included in the fit, which is for the case 
Mbh/M = 0.0075. The dotted line shows the profile without 
the contribution from the cusp, and in this figure the "cusp 
radius" is defined to be the radius at which the contributions 
from the core and cusp are equal. 
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computed density around the black hole increases with in- 
creasing black hole mass, this approach introduces a mass- 
dependent bias which leads to an underestimate of the 
true core radius for large black hole masses. Baumgardt 
et al. (2005) adopted a more observational procedure, and 
estimated the core radius as the radius at which the sur- 
face density drops to half its central value; but this is also 
biased in the same way. Therefore we have taken a dif- 
ferent approach, in which we fit the density distribution 
in the simulation by a template which describes the cusp 
around the black hole embedded in a core. Specifically, 
we fitted the density at radius r by 

p(rH(ar- 7 / 4 + 6)(l + r 2 /r c y, (7) 

where a,b,r c ,d are parameters. For a = (no cusp) this 
resembles closely the so-called EFF model for the distri- 
bution of density in a star cluster with a core (Elson et 
al., 1987; Mackey & Gilmore, 2003). The fit is not good 
at large radii; by experimenting, we found that the mean 
square residual was minimised if data outside about 1r% 
were rejected (Fig.l). 

Figure 2 shows fits of our theory to the data we have 
obtained in this way. The points on the left (crosses) come 
from one simulation, and plot the evolution of r c /rh as the 
black hole grows. The early points in particular do not fit 
the predicted power law very well, and we consider that 
these correspond to the period during which equilibrium 
has not yet been established between the flux of energy 
at the half-mass radius and that provided by the cusp. 
The points on the right (boxes) represent four different 
simulations, with different initial black hole mass, late in 
the evolution. The agreement with our theory is better, 
but the theory depends implicitly on a homology assump- 
tion, and the slightly discrepant slope may be due to small 
departures from homology; this is not unexpected, given 
that the values of r c /r^ extend up to about 0.8. 



Fig. 2. Dependence of core radius (in units of the half-mass 
radius) on black hole mass (in units of the total cluster mass). 
The line is the graph of r c /r h = 4.3(M bh /M) 3 / 4 (cf. eq.5), 
which is the best fit. Boxes: data from the four profiles in 
Baumgardt et al. (2004a); crosses: data at various times dur- 
ing the run with M(,h(0) = 266Mq (though the initial data 
point is not plotted). The error bars are the standard error 
provided by the fitting package, and include the contribution 
of sampling errors in the individual density measurements. 
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The important result to be obtained from this figure is 
that the core radius increases with black hole mass, which 
is the opposite of the conclusion that would be drawn 
from Table 1 of Baumgardt et al. (2005) 2 . There the core 
radius was estimated from the radius at which the surface 
density drops to half its central value (see above), and 
the fact that their paper considered clusters with a stellar 
mass function is not relevant in the present context. Trenti 
et al. (2007) found that the dependence of core radius on 
Mf,h is flatter than we predict, but there the quoted core 
radius was estimated as in iV-body models (see above). 

5. Discussion 

In this Letter we have argued that a balance is reached 
between the energy production of the cusp and the energy 
required for expansion of the entire cluster. Previously 
Marchant & Shapiro (1980) considered that there should 
be a balance between energy production of the cusp and 
that required for expansion of the core. This, however, 
leads to the relation r c /r cusp ~ constant, independent of 
the mass of the black hole and cluster. This is certainly 
incorrect when the the black hole mass is very small. 
Shapiro (1977), McMillan et al. (1981) and Duncan & 
Shapiro (1982) also, in effect, equated the generation of 
energy in the cusp to the energy required to expand the 
core; in fact, however, most of the energy passes through 
the core to the half-mass radius, as it provides the energy 
needed for the expansion of the entire cluster. Baumgardt 
et al. (2004a) equated the energy generation to that re- 
quired at rh, but assumed a constant value of r^/r c . Yuan 

2 We also note that the last value of Rc / Rh,pio m that Table 
should be 0.07. 
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& Zhong (1990) and Dokuchaev (1991) took a homology 
model for the stellar system, which was thus characterised 
by a single scale radius. 

The theory we have presented depends on a number of 
assumptions. 

1. It applies only to systems with stars of equal mass, 
though there is no reason to suppose that it does 
not extend, with suitable modification of detail, to 
systems with a realistic mass spectrum. 

2. It also depends, though in an inessential way, on 
the assumption that v 2 ~ v\ . In fact Y-body mod- 
els, albeit with a mass spectrum (Baumgardt et al., 
2005), show that the v 2 varies by only about 30% be- 
tween rh and O.lr/, . If instead we had assumed that 
v 2 oc r~^, where these results imply that (5 ~ 0.1, our 
result (eq.(5)) would have changed from a power law 
index of 3/4 to 6/(8 - 3/3) ~ 0.78. The velocity dis- 
persion profiles of some clusters require larger values 
of/3~0.3, however (e.g. u> Cen (Meylanet al., 1995), 
M15 (McNamara et al., 2003)). 

3. The theory applies only when sufficient time has 
elapsed for achieving a balance between the energy 
produced in the cusp and the flow of energy across 
the half-mass radius. In almost all galactic nuclei, 
and even in many globular clusters (Trenti, 2006), 
the relaxation time would be too long for the condi- 
tions assumed in this Letter to be established. 

4. The theory does not apply to systems where some 
other mechanism, e.g. interactions of primordial bi- 
naries, provide energy more efficiently. 

5. It does not apply to systems that are so dense (in- 
cluding many galactic nuclei) that the structure of 
the cusp is dominated by the role of physical col- 
lisions, and the power law is altered (Rauch, 1999; 
Duncan & Shapiro, 1983). 
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